Enhancement of pairing in a boson-fermion model for coupled ladders 
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Motivated by the presence of various charge inhomogeneities in strongly correlated systems of 
coupled ladders, a model of spatially separated bosonic and fermionic degrees of freedom is nu- 
merically studied. In this model, bosonic chains are connected to fermionic chains by two types of 
generalized Andreev couplings. It is shown that for both types of couplings the long-distance pairing 
correlations are enhanced. Near quarter filling, this effect is much larger for the splitting of a pair 
in electrons which go to the two neighboring fermionic chains than for a pair hopping process. It is 
argued that the pairing enhancement is a result of the nearest neighbor Coulomb repulsion which 
tunes the competition between pairing and charge ordering. 

PACS numbers: 74.20.-z, 74.45. +c, 74.81.-g, 71.10.Fd 



I. INTRODUCTION 

Charge inhomogeneities are an ubiquitous feature in 
strongly correlated electron systems. One of these in- 
homogeneities is the stripe phase present in some un- 
derdoped cupratesi In this phase, it has been recently 
suggested by theoretical— and experimental- studies that 
both stripes and intervening spin regions may be mod- 
eled as two-leg ladders. Charge inhomogeneity can also 
be originated by the structure of the materials. In the 
layered compound /3-Nao.33V2 0s, the V-0 planes con- 
sist of two-leg ladders separated by zig-zag chains^ This 
compound undergoes a transition from a charge-ordered 
stated to a superconducting state under pressured 
Within a purely electronic mechanism for superconduc- 
tivity in this material it is tempting to associate the for- 
mation of pairs to the ladder unitss£ Similarly, in the 
compound Sri4_ x Ca :r Cu2404i, the Cu-0 planes consist 
of coupled ladders forming a trellis lattice. This material 
becomes also superconductor under pressure^ and it also 
presents a competing charge ordered phased 

Then, it seems natural to study these inhomogeneous 
systems with Hubbard or t-J Hamiltonians defined on 
coupled ladders or quasi-one dimensional structures, with 
in general different fillings. After all, single ladders con- 
tains pairing, pseudogap, charge ordering (CDW)ii^ 
However, these models are quite difficult to analyze, both 
by analytical or numerical techniques. To make this 
problem more manageable, effective models of bosons 
and fermions can be derived and studied. These mod- 
els may be obtained after a basis change from the site 
to the dimer basisii*^ or to the plaquette basisi^ and 
eventually projecting out some states of the new basis. 

Simplified models where the hard-core bosons describ- 
ing triplet excitations have been eliminated have also 
been studied recently^ These models were originally 
proposed to describe bipolaronic superconductivity^ and 
later used to study the pseudogap phase in the cuprates 16 
where typically bosons represent preformed pairs. In the 
present work, as a difference with those earlier studies, a 
model in which bosonic and fermionic degrees of freedom 



are spatially separated is considered. 17 

In the case of the stripe order, as suggested by one of 
the main theories about this phased preformed pairs on 
the stripes would be described by hard-core bosons while 
unpaired charge in the intervening regions would be de- 
scribed by fermions. Although pairing on stripes is a con- 
troversial issued it is still interesting to examine a pos- 
sible enhancement of superconductivity due to proximity 
effect within a simple lattice model. In the same way, the 
compounds with structural ladders could also be mod- 
eled by spatially separated bosons and fermions. Bosons 
would describe pairs on ladders and fermions unpaired 
electrons on the zig-zag chains. In both cases, a rela- 
tively simple boson-fermion model could give qualitative 
features about, for example, the competition between 
superconductivity and CDW, the effect of strong elec- 
tron correlations, and the effect of applied pressure which 
leads to the modification of couplings, site energies^ 

In general, these effective models will contain fermion 
and boson hopping terms together with some additional 
terms mixing bosons and fermions. One of the most im- 
portant and interesting mixing term is a generalized An- 
dreev coupling which describes the breaking of a pair in 
two electrons or the reverse process. In Section [0] the 
model studied in this paper, formulated on a system of 
alternating bosonic and fermionic chains, is derived from 
microscopic models using a projection technique. Results 
obtained by exact diagonalization are shown in Section 
IIIII Finally, the relevance of these results to various phys- 
ical systems is discussed in Section llVl 



II. EFFECTIVE MODEL 



The specific model studied in this paper is formulated 
on a system of alternating bosonic and fermionic chains 
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FIG. 1: (a) Inhomogeneous system of coupled chains with 
bosonic (thick lines) or fermionic degrees of freedom (thin 
lines), (b) Two different processes of creation/annihilation of 
pairs: splitting (top) and pair hopping (bottom), (c) deriva- 
tion of an effective model by changing to a dimer basis fol- 
lowed by projection, (d) effective couplings (defined in the 
text) as a function of the inter-dimer hopping, J = 0.4, Vb = 
and t = 1 in the original 12-site cluster, (e) Same as (d) for 
U = -8, Vb = 2 and t = 1. 



as shown in Fig.^a). The Hamiltonian is defined as: 
H = - t b J2 (4bj + h.c.) ~tf J2 (4 c i- + h - c -) 

<i,j> <i,j>,a 

+ U/]n^nji + y ^ eimejuj 

i <i,j> 

+ ^eiUi + (i>i CjtCfej, +h.c.) (1) 

where, cj CT creates an electron with spin a at site j, 
rij a = Cj a Cj a ; rij — rij^ +n>ji for fermions or nj = b^bj for 
hard-core bosons. The charge is ej = 1 (2) for fermions 
(bosons), tb and tf are the hopping integrals along the 
chains for bosons and fermions respectively; U is the on- 
site Coulomb repulsion on the fermion sites; V is the 
Coulomb repulsion on nearest neighbor (NN) sites, act- 
ing between paired and/or unpaired electrons; 6j is the 
onsite energy at site i which we take = e on the bosonic 
chains, = on the fermionic chains. Thus, e contains 
the binding energy. 

The last term of the boson- fermion Hamiltonian Eq. 
is a generalized Andreev coupling, transversal to the 
chains, where i is a site on a bosonic chain and j, k 
are the NN sites of i on the fermionic chains. Two differ- 
ent processes will be studied (Fig. ^b)). The one at the 
top corresponds to a splitting of a pair into two electrons 
located in the two fermionic neighboring chains (A^s) 
while the one at the bottom corresponds just to a pair 
hopping to one of the two neighboring fermionic chains 



In order to understand the physical origin of these two 
types of Andreev coupling let us consider an extended t-J 
or Hubbard model on the 12-site system of Fig. [He) top 
with 2 electrons. By performing a site to dimer change 
of basis, and projecting out the one-electron (fermionic) 
dimer states on the two inner dimers and the double oc- 
cupied (bosonic) states on the four outer dimers, one gets 
the effective 6-site boson-fermion model of Fig.^c) bot- 
tom. The effective Hamiltonian is given by the standard 
formula: 

H eff =VH V- VH Q QHq q_ Eq 2 W (2) 

where V is the projection operator on the subspace of 
retained states, Q is the projection operator on the sub- 
space of the eliminated states, TCo is the original Hamil- 
tonian, Ho^o = Eq^q, and H ef fV^ = E V^ M Vari- 
ants of this procedure were repeatedly performedii'i2ii& 
to obtain effective models from the Hubbard or t-J mod- 
els, specially retaining triplet states and projecting out 
fermionic states. Similar studies but retaining fermionic 
states, although using a different projection procedure^, 
have concluded that the most important interactions not 
involving triplets are the ones contained in Hamiltonian 
(JTJ. Although Eq. (0) is usually analytically calculated 
using second order perturbation theory, for the 12-site 
cluster of Fig. ^c) it could easily be numerically solved 
using standard matrix inversion subroutines^ In the 
mapping shown in Fig.^c), the projection step also leads 
to second neighbor fermion-fermion and three-site inter- 
actions in the horizontal direction which are also negligi- 
ble in first approximation. 

Let us consider first that H.q is an extended t-J model. 
Let us assume that we have the usual t, J couplings and 
a NN Coulomb repulsion Vq on the ladders, and ti n t er , 
Jinter, V lnte r between ladders, with J in ter/J = {Unter/t) 2 
and Vi n t e r/Vo = Unter/t- It is reasonable to assume that 
the main effect of applying pressure is to modify the in- 
ter ladder interactions, which are here related to U n t e r- 
The values of the effective couplings tb, tf, Xa, s and 
are shown in Fig.^d) for J/t = 0.4, VoJ/t — 0, as a func- 
tion of ti n ter/t. The most important feature is that A^ ;S 
is clearly larger than \A,hop- It should be noticed that 
the sign of Xa, s corresponding to the process in Fig.^b) 
is opposite to the process in which the up spin electron 
goes to the left chain and the down spin electron goes to 
the right chain. This is related to the fact that the elec- 
tron pair on a ladder rung form a singlet. One should 
also notice that tb and tf take values close to each other. 

On the other hand, let us assume that we have an ex- 
tended Hubbard model with Uo < 0, Vq > 0, on the 12- 
site cluster of Fig.njc). This attractive Uq could describe 
a phonon mediated pairing. In this case, one obtains the 
effective parameters shown in Fig. ^e). It is now ap- 
parent that the pair hopping type of Andreev coupling 
dominates over the pair splitting type. It should also be 
noticed that in this case tb is much smaller than tt and 
hence it could be neglected which is precisely what has 



been done in the earlier literatureiai^ although it was 
included in the studies related to the pseudogap phase 
in cuprates»i& It is possible to think then that the case 
considered in Fig^d), with a dominance of \a,s, cor- 
responds to a strongly correlated electron physics with 
a likely d-wave pairing, while the situation of Fig^c), 
would correspond to a more conventional, s-wave, type 
of superconductivity. 

In order to support this interpretation, the probabil- 
ity of double occupancy and the probability of having a 
singlet on a rung, properly normalized (i.e., the sum of 
all the possible configurations on a rung equal to 1) was 
computed. In the case of the extended t — J model, the 
probability of electrons forming a rung singlet is much 
larger than the probability of electrons going to double- 
occupied sites. The reverse situation occurs for the at- 
tractive Hubbard model. It is instructive to consider also 
a model which interpolates between the Hubbard and 
the t-J model, the so-called t-J-U model, obtained from 
the t-J model by relaxing the no double-occupancy con- 
straint but including an onsite Hubbard repulsion. For 
an intermediate situation, for example, J = 1, t = 1, 
U = 0.5, V = 0, there is a crossover from the splitting 
to the pair-hopping types of Andreev coupling as ti nter 
is increased consistent with a crossing from singlet to 
double-occupancy order. That is, applied pressure can 
change a s-wave pairing into a d-wave pairing. 

The proposed model is far more general than the 
simple "derivation" schematically shown in Fig. QJc). In 
the first place, a boson does not necessarily represent a 
pair on a ladder rung. In fact, it has been shown that 
on ladders, pairs are located along plaquette diagonals 
or on more distant sites depending on coupling values^ 
In general, this boson-fermion model is applicable to any 
compound containing quasi-lD units bearing some kind 
of pairing. The site energy ej would be in general deter- 
mined by the binding energy of electrons on these quasi- 
lD units as well as by a potential coming from the whole 
structure of the material, which can be modified by ex- 
ternal applied pressure or by internal chemical substitu- 
tion, as in Sri4_ x Ca x Cu2404i, where Ca doping leads 
to transfer of holes from the chain to the ladder planes. 
The boson-fermion model could provide insights to pre- 
dict the effects produced by these kinds of perturbations 
on a given compound. 



III. NUMERICAL RESULTS 

A. 3 x L cluster, quarter filling 

Model Eq.JTJ was studied by exact diagonalization 
(Lanczos algorithm) on 3xi (L = 6, 8) clusters with peri- 
odic (open) boundary conditions (BC) along (across) the 
L-site chains. The central L-site chain has bosonic oper- 
ators while the two external chains contain the fermionic 
ones, tf is adopted as the unit of energy. In the above 
mentioned basis change, the effective parameters U and 




FIG. 2: Results for the 3 x 8 cluster, U = V = 2, e = 0, 
n = 0.5, as a function of if, and \a- Boson density in the 
central chain: (a) pair splitting, (b) pair hopping. Pair-pair 
correlations at the maximum distance along the central chain: 
(c) pair splitting, (d) pair hopping. 



V result roughly half the NN Coulomb repulsion of the 
original model (Vq), both for the limits of infinite and 
zero values of the Hubbard on-site repulsion of the orig- 
inal model (Uq). Hence U = V in the following. Various 
properties, specially those related to superconductivity, 
were computed as a function of tb, e and Xa- 

Figure |21 shows results for the 3x8 cluster, 6 = 0,11 = 

V = 2, at quarter filling (n = 0.5). The first feature to 
notice is that the relative occupation of the fermionic and 
bosonic chains depends on the parameters of the model. 
In Figs. EI a) and (b), the boson density in the central 
chain, 6b = (rib)/L (L = 8 in this case) is shown for Xa, s 
and \A,h respectively. In both cases, 5b increases as tb is 
increased. This is expected since electrons would move to 
the central chain to gain kinetic energy. For tb = 1 there 
is a level crossing with a sudden increase of 8b- On the 
other hand, 8b slowly decreases with \A,h in the interval 
shown, while its behavior with \a, s is non monotonic. 

The central quantity of the present study is the boson 
correlation at the maximum distance, P Tmax = (b r bo) 
along the central chain. This correlation, which in the 
current model has the meaning of pairing correlation, 
is a measure of quasi-long range superconducting order 
on the bosonic chain. The results for P rmax are shown 
in Fig. 12c), (d) for A,4, s and \A,h respectively. In both 
cases, for a fixed value of tb, Pr max shows an enhance- 
ment as the Andreev coupling increases. The curves of 
P Tmax vs. \a are shifted upward as tb increases, as ex- 
pected. The most important feature of these results is 
that the enhancement with A^ is much stronger for the 
case of pair splitting than for the case of pair hopping. 
Notice also that for pair splitting, for if, fixed, the behav- 
ior of Pr max and the one of 8b are unrelated. It should be 
emphasized that in the region where P Tmax is enhanced, 
the pairing correlations as a function of distance, (fej&o), 
have a monotonically decreasing behavior corresponding 
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FIG. 3: Results for the 3 x 8 cluster, U = V = 2, e = 0, n = 
0.5, as a function of tt and \a- Maximum charge structure 
factor: (a) pair splitting, (b) pair hopping. Current-current 
correlations at the maximum distance along the fermionic 
chains: (c) pair splitting, (d) pair hopping. 



to true long-distance pairing. A non-monotonic behavior 
would be indicative of phase separation or CDW. 

For U = V = 0, energy, boson occupancy and pair- 
ing correlations are identical for Xa, s and Xa,h- This 
comes from the fact that the respective Hamiltonians 
are related by the transformation — T~ 1 'H^T, 

where T relates the two processes depicted in Fig. ^b) 
as Hh = T~ 1 H S T, and T _1 = T T . Quantities defined 
solely in terms of bosonic operators are preserved by this 
transformation. The important result is that in this case, 
Pr max is not enhanced by Aa, although it is considerably 
larger than for U = V = 2. 

In order to characterize the physics of this model more 
completely, the static charge structure factor C(q) and 
the current-current correlations at the maximum distance 
along one of the fermionic chains were computed. The 
first quantity is indicative of CDW while the second one 
is related to the conduction of the fermionic chains^ For 
the same parameters as before (U = V = 2, e = 0, n — 
0.5), C(q) presents a maximum at q max = (37r/4, 2tt/3) 
in the whole range of % and Xa examined. Figures 
03a), (b) show C m ax = C(q mal ) for A^.s and X A: h respec- 
tively. It can be seen that this quantity is suppressed, 
particularly by Xa.s- The charge structure factor for the 
whole cluster behaves in a similar way than the one com- 
puted from the charge-charge correlations along a single 
fermionic chain. Now, for U = V = 0, C(q), as it was 
found for P rmam , is roughly independent of Xa- The same 
behavior is also found for other clusters and densities con- 
sidered below. It is possible then to sum up these fea- 
tures by suggesting that Xa.s works against the tendency 
to CDW, favored by V, thus leading to an enhancement 
of long distance pairing. 

The current operator between sites i and i+y is defined 

as usual as jy,i = ^J2a( c i+ycr Cicr ~ h.c), and current- 
current correlations at the maximum distance as ay = 
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FIG. 4: Results for the 3 x 6 cluster, t b = 0.8, n = 0.67 and 
0.89. (a) Boson density and (b) pair-pair correlations at the 
maximum distance along the central chain for U = V = 1, 
e = —0.5, as a function of \a- Results for pair splitting 
(hopping) are indicated with circles (triangles). In (c) and 
(d), the same quantities as in (a) and (b) respectively are 
shown for U = V = 0.5, e = -2. 



(jy,r ma:c jy,o)- Results for oy max along a fermionic chain 
are shown in Fig.|3fc),(d) for Xa.s and Aa,/i respectively. 
It can be seen that o> max is suppressed in both cases 
although this effect is much larger for Xa, s than for Xa,H- 
This behavior indicates that the effect of A^ is to favour 
the conduction mainly through the bosonic chain. 



B. Other clusters and fillings 

It is also important to determine if the behavior shown 
in Figs. |5] and is also present at larger electron fill- 
ings, specially because some possible applications of the 
present model, for example Sri4_ 2; Ca a ;Cu2404i, corre- 
spond to systems close to half-filling. As electron fill- 
ing increases from n = 0.5, the dimension of the Hilbert 
space increases rapidly and it soon makes this problem 
very hard for exact diagonalization. Hence we have to 
limit the study to the smaller 3x6 cluster. 

Results for the 3x6 cluster with 12 electrons (n — 0.67) 
and 16 electrons (n — 0.89) are shown in Fig. 0] For 
U = V = 1, e = —0.5, tb = 0.8, the boson density varies 
slowly with X A (Fig. Ufa)). At n = 0.67, S b - 0.35, 
which implies an identical charge density on fermionic 
and bosonic chains. At n = 0.89, 5b ~ 0.5 imply- 
ing a larger charge density on the bosonic than on the 
fermionic chains. The overall behavior of the pairing 
correlation at the maximum distance along the central 
chain, shown in Fig. Ilfb), is the same as in Fig. i.e., 
there is an enhancement of P Tma7 . with A^ which is more 
pronounced for pair splitting process. For these values of 
the parameters, C max , peaked at q max = (tt, tt) is mono- 
tonically suppressed by A^. 

To obtain larger charge density on the bosonic chain, 
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FIG. 5: Pair-pair correlations at the maximum distance along 
the central chain versus the absolute value of the energy of 
the Andreev term. Results for the 3x8 cluster, U = V — 2, 
n = 0.5, tb = 0.6 (a) e = 0, (b) e = -0.5. Results for the 
3x6 cluster, U = V = 1, tb = 0.8, e = -0.5, (c) n = 0.67, 
(d) n = 0.89. 



the couplings U — V — 0.5 and e = —2 were studied; 
tb = 0.8 as before. It can be seen in Fig. 0Jc), that 
at n = 0.67, Sb becomes slightly larger than 0.5 and at 
n = 0.89, Sb ~ 0.6. For both global fillings, the charge 
density on the bosonic chain is approximately twice the 
one on the fermionic chains. Notice that now P Tmax 
(Fig. |2[d)) is larger than the one shown in Fig. 0{b). 
This larger value, for n — 0.89, could be related to the 
behavior of simple hard-core boson (or spinless fcrmion) 
chain with NN repulsion, where superconductivity is sup- 
pressed at half-filling (5b = 0.5), which is the case for 
U = V = 1 and e = -0.5 (Fig.gfa)). On the other hand, 
the opposite happens for the case of n = 0.67, indicating 
that the Andreev coupling changes the physics of an iso- 
lated bosonic chain. However, at A^i < 0.2, the pairing 
correlations as a function of distance has a non mono- 
tonic behavior, signaling CDW. It should also be noticed 
that for all cases in Fig. 0fd), there is a saturation and 
further decreasing of -Pr max for larger A^- This may in- 
dicate that the behavior shown in Figs.[2Ic) and^Jb) is 
mostly a property of low bosonic density (Sb < 0.5). It 
should be stressed that by going from e = —0.5 to —2, 
with U = V = 1 fixed, the changes are smoother than 
by going from U — V — 1 to U — V — 0.5, keeping 
e = —2 fixed. In this parameter space, results interpo- 
late smoothly between those of Fig.|3Ja),(b) and those of 
Fig. He), (d). 

It is tempting to relate the relatively small effect of pair 
hopping on pairing correlations to its possible suppres- 
sion by the onsite Coulomb repulsion on the fermionic 
chains. This is actually not the case. One should notice 
first that pair splitting is also affected by such repulsion 
since an electron could be already present in one or both 
of the final sites on the fermionic chains. For noninteract- 
ing electrons simple combinatorics lead to the result that 



double-occupancy is more likely on the splitting than on 
the pair processes for electron densities larger than ~ 0.6. 
A similar effect caused by the NN repulsion V is more 
difficult to predict. Alternatively, it is possible to com- 
pute the contribution of the Andreev term to the total 
energy, Ea, as a measure of how much this term is actu- 
ally "working" . It may be convenient then to plot P Tma] . 
as a function of Ea, rather than as a function of the bare 
parameter Xa- This is done in Fig. for some typical 
cases of Figs. [21 and 0] It can be seen that the qualita- 
tive behavior of these figures is not modified. Only in 
Fig-E^d), corresponding to a density x — 0.89, there is a 
jump in P rmax for the pair hopping case but this occurs 
at a rather large value of the bare coupling, Xam = 1-2. 
Figs. 0a) and (b) allows the comparison between two dif- 
ferent values of e. For both types of Andreev couplings, a 
smaller (negative) value of e gives a smaller enhancement 
of Pr max , an even a suppression in the case of the pair 
type of Andreev coupling. 

Computer limitations make it difficult to go to larger 
clusters in order to assess finite size effects but it is 
possible to study clusters with different geometry. The 
2x12 cluster was considered to estimate finite size effects 
on the pair hopping between the fermionic and bosonic 
chains. Results for the same parameters as in Fig. [21 
show that P m ax takes values very close to those found 
for the 3x8 cluster with a very similar (small) depen- 
dence with A^/j. This asymmetric two-chain system is 
essentially the same studied by Le Hur— by bosoniza- 
tion techniques, although in this work the boson medi- 
ated pairing of unpaired electrons is mainly investigated. 
Finally, the cluster with 4 coupled chains of length 6, 
with periodic BC also in the transversal direction, with 
12 and 14 electrons, and the same couplings as in Figs. [21 
and [3 was considered. The overall behavior is the same 
as that depicted in those Figures -in particular P Tma7 . is 
much more robust for Aa, s than for \a,k- with the addi- 
tional feature of an enhancement of pairing correlations 
also in the direction perpendicular to the chains^ For 
U = V = 0, as for the 3 x L clusters studied above, 
C(q) and longitudinal P Tmaus are almost independent of 
Xa although transversal pairing correlations are trivially 
enhanced by A^. 



C. A "toy" model of stripes 

With respect to the application of the present model 
to the stripe phase in the cuprates, one should take into 
account that the t-J model considered as a microscopic 
model from which the couplings in Fig. ^d) were de- 
rived, was meant to be the strong-coupling limit of the 
one-band Hubbard model. In these models pairs involve 
electrons with opposite spins. In the t-J model, as ap- 
plied to cuprates, pairing involve doped holes which are 
described by singlets^ Although a derivation of an effec- 
tive model from this version of the t-J model is possible, 
it is instructive to use the already obtained results to 
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FIG. 6: Results for the 4 x 6 cluster, n = 0.25, U = V = 0, 
tb = tf = 1, Xa,s = 0.33, \a,Ii = 0.06. (a) Boson density, and 
(b) pair-pair correlations at the maximum distance along a 
bosonic chain (filled triangles) and perpendicular to the chains 
(open circles) versus the absolute value bosonic site energy. 



study a "toy" model of stripes. To do this, a fermion 
should have the meaning of a doped hole, and the half- 
filled state of the cuprates should be the "vacuum" of 
model |QJ. From Fig. IHd), the approximate values of 
the couplings are tb = tf = 1, Xa.s = 0.33, \A,h = 0.06, 
U = V = 0. On the 4x6 cluster described above, in 
order that the stripes be at a linear filling of one quarter, 
there should be 6 doped holes, corresponding to a doping 
on the original cluster of 0.125. In this study the vari- 
able is the site energy at the stripes, e, which may also 
depend of various mechanisms such as structural details, 
phonons^S Results are shown in Fig. |SJ The main con- 
clusion is that pairing is enhanced in both longitudinal 
and transversal directions, even though doped holes are 
increasingly localized at the "stripes" . 



pounds in which pairing takes place in quasi-lD struc- 
tures such as two-leg ladders. Starting from a micro- 
scopic model, an exact projection procedure on a small 
cluster suggests that a pair splitting kind of Andreev 
process is related to the physics of repulsive U systems, 
characterized by d-wave pairing, while a pair hopping 
Andreev process is more related to a negative U kind 
of physics leading to s-wave pairing. This elementary 
projection also gives indication of how the effective cou- 
plings are changed with pressure. The values of these 
couplings for a specific compound should be obtained 
from a realistic, in general complex, microscopic model, 
and in this case an exact diagonalization procedure deal- 
ing with much larger clusters than the ones of Fig. Qfc) 
should be used. Although in principle both types of An- 
dreev couplings are going to be present in the effective 
model irrespective of the nature of the pairing, the pur- 
pose of the present study was to determine the more gen- 
eral and important properties of those processes taken 
separately. The conclusion was that, close to quarter- 
filling, a pair splitting process is more efficient to en- 
hance long-distance pairing than pair hopping from the 
superconducting to the non superconducting chains. So 
far this result would suggest that if the effect of pres- 
sure translates in an increasing of A,4 in, for example, (3- 
Nao.33V2C>5, then the presence of pairing would be more 
likely the result of strongly correlated electron physics in 
these compounds. More detailed predictions would re- 
quire to determine the effective couplings more precisely 
as discussed above. In any case, even at this general level, 
more predictions for example for ARPES experiments 
could be obtained by studying dynamical properties^ 



IV. CONCLUSIONS 

In summary, a model of coupled bosonic and fermionic 
chains was proposed to describe the physics of com- 
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